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Abstract 

We present a zero decomposition theorem and an algorithm based on Wu's 
method, which computes a zero decomposition with multiplicity for a given zero- 
dimensional polynomial system. If the system satisfies some condition, the zero 
decomposition is of triangular form. 
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1 Introduction 

Since Wu's method [8 1 was invented, many algorithms for zero decomposition of poly- 
nomial systems have been proposed fl] H] [6] [7] [9] . The so-called zero decomposition 
is to decompose the zero set defined by a polynomial set into a union of finite many 
constructible sets. Usually, these constructible sets are defined by some triangular 
sets, i.e., the equations defining each constructible set are in triangular form, e.g., 
fi(xi), f2(xi,X2), ■ ■ ■ , f n {x\, ...,x n ). The zeros of such equations can be described 
and solved easily, so the methods of zero decomposition have been applied success- 
fully to many problems with various backgrounds, especially those in geometry. 

All the existing methods for zero decomposition focus on the decomposition of the 
set of zeros and do not care about the multiplicities of (isolated) zeros. For a zero- 
dimensional polynomial system, we are interested in getting its zero decomposition 
with multiplicity. More precisely, for a zero-dimensional polynomial system S, we 
hope we can compute finite many pairs [Tj, gi] such that 

MZero(S') = (J MZero(iy 5i ), 

i 

where MZero(-) stands for the set of zeros (counted with multiplicity) and MZeio(Ti/ gi) 
{£ : V/ e Tj, /(£) = 0, 0} (counted with multiplicity). Herein, Ti is a poly- 

nomial set, gi is a polynomial and for all i ^ j, MZero(T,/ gi) n MZero(Tj / gj) = 0. 
That is to say, for any zero p of S, there exists only one pair [Tj,gj] such that p is a 
zero of Tj and gj (p) ^ and the multiplicity of p as a zero of S is the same as that of 
it as a zero of Tj. We will call such decomposition keep multiplicity. 

For a known zero p of a zero-dimensional system S, one can compute the multiplic- 
ity of p through Grobner bases computation [2|. Furthermore, one can even compute 

'Corresponding author, email: xbc@math.pku.edu.cn 



1 



the multiplicity structure of p by the method in [ 3 1 using dual space. Obviously, our 
interest is different from those in |2| and |3|. We hope we can compute a zero decom- 
position which keeps multiplicity for a zero-dimensional system by modifying Wu's 
method. In that way, we will not add too much additional heavy computation to the 
original method. Our interest is also different from that in [5 1, where a new definition of 
multiplicity based on non-standard analysis was given and was proven to be equivalent 
to the classical definition. Zero decomposition is not a topic in J5). 

We hope that a zero decomposition keeping multiplicity is also triangular. If so, 
the multiplicities of zeros of each component can be read out directly ifTUl . In this 
paper, we give a zero decomposition theorem and an algorithm based on Wu's method, 
which computes a zero decomposition with multiplicity for a given zero-dimensional 
polynomial system. If the system satisfies some assumption, the zero decomposition is 
of triangular form. 

The organization of the paper is as follows. In Section 2, some preliminary knowl- 
edge about the definition of multiplicity is reviewed. Section 3 proves a theorem of 
zero decomposition with multiplicity for given zero-dimensional polynomial systems. 
Section 4 presents the algorithm of zero decomposition with multiplicity and points out 
the condition for the algorithm to output triangular decomposition. We discuss some 
possible improvements on the algorithm in Section 5. 

2 Preliminaries 

First, let's recall the definition of local (intersection) multiplicity. We follow the nota- 
tions in |2|. Although some notations and definitions can be stated in a more general 
way, we restrict ourselves to the ring C[X] = C[x\, . . . , x n ] since we are interested in 
the complex or real zeros of zero-dimensional polynomial systems. 

For p = (771, ... , r) n ) £ C™, we denote by M p the maximal ideal generated by 
{x\ — rji, . . . ,x n — rj n } in C[X], and write 



It is well-known that C[X]m p is the so-called local ring. 

Definition 1 Let I be a zero-dimensional ideal in C[X], and assume that p G Zero(Z), 
the zero set of I in C. Then the multiplicity of p as a point in Zero(Z) is defined to be 



the dimension of the quotient space C[X]m p /IC[X]m p as a vector space in C. 

There is another equivalent definition of local multiplicity using dual space. 
For every index array j = . ■ ■ , j n ] £ N n , define the differential operator 




dim c C[X] Mp /IC[X} Mi 



9j = a 



= d 



I hin 
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X' 



•••X' 



n 



and the order of dj is \dj \ = Yl 3i- 

1=1 
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Consider a system of m polynomials, {fi(X), . . . , f m (X)}, in n variables with an 
isolated zero £ 6 C", where m > n. Let / = (fi(X), . . . , f m (X)). A functional is 
defined at £ 6 C" as follows: 

dj[$ : C[X] -> c, 

where 

dM(f) = (djfm, V/6C[I]. 
All functionals at £ that vanish on / form a vector space D^(I), 

D ( (I) ^{v=J2 v J d AC] ■■ »(/) = 0,V/ e 1} 

where Vj £ C. The vector space D^(I) is called the c/wa/ space of / at £. For a — 
0, 1, . . . , D^(I) consists of functionals in D^(I) with order bounded by a. 

Definition 2 [3] The local multiplicity of zero £ of a zero -dimensional ideal I C C[X] 
is m if the dual space D^(I) is of dimension m. 

Remark 1 According to [3$, Definition\2\is equivalent to Definition^ 

Proposition 1 Q Let a be the smallest a such that dim(-D|") = dini(Z)^ +1 (/)), then 
D^{I) = Dg. Furthermore we have 

a <m = dim(D^(/)). 

Lemma 1 Assume £ is an isolated zero of I and the local multiplicity of £ is m. If 
£ S Zero(g<(X)) where g(X) e C[X], then for any v G D^(I), any h(X) E C[X] and 
l>m,we have v{h{X)g{X) 1 ) = 0. 

Proof. If dj = d n ... jn e D ( (I) = D%(I), \dj\ < a < mby Proposition Q] Be- 
cause g(X) l -\ d i I is a factor of dj{hg l ),dj[i](h(X)g(X) 1 ) = since £ e Zem(g{X)). 
Because any v <E D^(I) is a linear combination of dj G D^(I) = D^(I), we are done. 

3 Zero decomposition with multiplicity 

Existing algorithms for triangular decomposition of zero-dimensional polynomial sys- 
tems do not care about the multiplicities of the systems' zeros. For example, Wu's 
method does not keep the local multiplicities of zeros. We will prove in this section 
that for a given zero-dimensional polynomial system S, there is an algorithm which 
computes finite many pairs [Ti,gi] where Tj is a characteristic set and gi is a polyno- 
mial such that 

MZero(S) = (J MZero(T i / 3i ) 

i 

as long as S satisfies some special conditions. Herein, MZero(-) is the set of zeros 
counted with multiplicity. Also, we have 

MZero(7y 3l ) n MZero(T.,- / gj) = if i ^ j. 
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We mainly follow Wu's method and make necessary modifications for our purpose. 
Throughout this paper, we use without definition some basic concepts and notations 
from [8 1 with which we assume the readers are familiar. 

For two polynomials f(x) and g(x), prem(/, g, x) is the pseudo remainder of / 
pseudo-divided by g with respect to (w.r.t.) x. For a polynomial f(X) and a tri- 
angular set T, prem(/, T) is the (successive) pseudo remainder of / w.r.t. T. If 
r = prem(/, T), it is clear that r G {/, T), the ideal generated by {/} U T. 

According to Wu's method, for a given polynomial system S, one can compute a 
so-called characteristic set of S, say C = [Ci(X), . . . , C r (X)\, such that 

r 

Zero(S') = Zero(C/ J) U (J Zero(S, J,-), (1) 

i=i 

where I L is the initial (leading coefficient) of Cj(X) for i — 1, . . . , r, J = Y[j=i Ij> 
Zero(C/J) = Zero(C) \Zero(J) and Zero (S, U) = Zevo(S U {/»}). Throughout this 
paper, we name the algorithm in [8 1 for computing the characteristic set C in (fl} as 

WuCharSet. 

The following lemma shows that the local multiplicity of each point in Zero(C/ J) 
is unchanged. That is to say, the first characteristic set computed by WuCharSet 
keeps multiplicity unchanged. 

Lemma 2 Assume S = {/i(X), . . . , f s (X)}, I = (S) is a zero-dimensional ideal, 
and C = [C\(X), . . . , C r (X)\ is a characteristic set of S computed by Wu's method 
such that (Q holds. Then for each p G Zero(C/J), the local multiplicity of p as a point 
in Zero(/) and the local multiplicity of p as a point in Zero ((C)) are the same. 

Proof. Assumep = (p±, . . . ,p n ) g Zero(C/J), andletM p = (x±— pi, . . . ,x n — 
p n ). According to Definition Q] the local multiplicity of p as a point in Zero(J) is 
dime C[a;i, . . . , x„\m p /IC\x\, . . . , x n ]M p and the local multiplicity of p as a point in 
Zero((C)) is dim c C[xi, . . . , x n } Mp /(C)C[xi, . . . ,x n ] Mp - It suffices to prove 

IC[xi, . . .,x n ] Mp = (C)C[xi,. . .,x n ] Mp . 

On one hand, (C)C[xi, . . . , x u ]m p Q IC[x\, . . . , x„]m p since (C) C /. On the other 
hand, prem(/i(Jf), C) = for each fi(X) G S since C is a characteristic set of S. 
Therefore, there exist oi, . . . , a r , qi(X), . . . , q r (X) such that 

f[l j ; i fi=j2q j (X)C j (X). 

Because p G Zero(C/J), Ij(p) ^ for all 1 < j < r. Thus, each Ij(X) is 
invertible in C[xi, . . . , x n ]M p - As a result, fj(X) G (C)C[xi, . . . , x n ]M p . Thus, 

IC[xi, . . .,x n ] Mp C (C)C[xi,. . .,x n ] Mp . 

That completes the proof. 

For those points in Zeio(S, Ii), their local multiplicities as points of Zeio(S,Ii) 
may be not equal to their local multiplicities as points of Zero(S'). Intuitively, if p G 
Zero(5, Ii) and its local multiplicity as a point of Zero (5) is mi, its local multiplicity 
as a point of Zero(S, I" 1 ' 2 ) should also be mi provided that > mi. And this is 
indeed true as implied by the following lemma. 
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Lemma 3 Suppose^ £ Zero(S', g{X)) is an isolated zero of S — {fi(X), . . . , f s (X)} 
and m is the local multiplicity oft; as a point ofZ£io(S). Then the local multiplicity of 
£ as a point ofZeio(S, g(X) 1 ) is m if I > m. 

Proof. Let / = (S) and V = (S,g(X) 1 }. According to Definition |2] it is 
sufficient to prove that 

Dz(I) = Ds(l'). 

Obviously, D^(I) D D^(I') according to the definition of dual space since / C I'. 

On the other hand, for any v € D £ (J) and / = £ d j (X)f j (X) + h(X)g(X) 1 E 

3=1 

I', 

S 

v(f(X)) = v(Y / d 3 (X)f ] (X))+v(h(X)g(Xy)=v(h(X)g(X) 1 ) 

3 = 1 

s 

because dj(X)fj(X) e /. According to Lemma[T] v{h(X)g(X) 1 ) = and thus 

3 = 1 

v e D^(I'). That completes the proof. 

By Lemma|2] Lemma|3]and eq.([T]i, the following theorem is obvious. 

Theorem 1 (Zero Decomposition with Multiplicity) Suppose S is a finite non-empty 
set of non-zero polynomials in C[X] and (S) is a zero-dimensional ideal. Assume that 
C = [Ci(X), . . . , C n {X)\ is a characteristic set of S, Ii is the initial of d(X) and 

i 

J{i) = Y[ Ik far 1 < i < n. There exists € Nfor 2 < i < n such that 

k=l 

n 

MZero(5) = MZero(C*/J(n)) U (J MZero(5 l / J(i - 1)), (2) 

i=2 

where Si = S U {1^'}, is a disjoint decomposition o/MZero(5). 

The existence of such e^s is clear. For example, we can let e, be the Bezout bound, 
or set e, > m where m — dime C[xi, . . . , x n ]/ (S) is the number (counted with 
multiplicity) of points in MZero(S'), which can be obtained by computing the Grobner 
basis of (S) |2). 

Additional, the characteristic set C in (2) satisfies (S) = (S U C), so we replace S 
by S U C sometimes to make sure our algorithm terminates. 

4 Algorithm and examples 

For a given S in Theorem[T] the characteristic set in eq.© can be computed through 
WuCharSet. Naturally, one may hope to apply WuCharSet to each Si recursively 
to obtain a triangular decomposition of MZero(5). However, if we directly apply 
WuCharSet to Si or Si U C recursively, the process may not terminate. 

We take the following strategy. Let S[ = S U C U {r;} where = prem(/ t ei , C). 
We replace Si by S[ in the recursive process. Because is reduced w.r.t. C, the rank 
of the characteristic set of S[ must be less than that of C provided that r,; ^ 0. 

However, we have to prove that this modification does not change the local multi- 
plicity of each point in Zero(Si/ J(i — 1)) = Zeio(S' i / J(i — 1)). 
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Proposition 2 Let the notations be as above. Then for each i > 2 and a point p G 
Zero(Si/J(i — 1)), the local multiplicity of p as a point of Xexo(Si) is equal to its 
local multiplicity as a point in Zero(S','). 

Proof. Assume that p = (pi, . . . ,p n ) and M p — (x\ — p\, . . . , x n — p n ). It is 
sufficient to prove that 

(S t )C[X] Mp = (S>)C[X] Mp . 

Since C is a characteristic set of S, (C) C (S). U is a polynomial in x\, . . . , x.^i be- 
cause (S) is zero-dimensional. Then there exist ni, . . . ,rti-i G Nand(7i(X), . . . G 
C[X] such that 

3=1 3=1 

On one hand, it is obvious that i\ G (Si), so 

(S t )C[X}M p 3 (^)C[X] Mp . 

On the other hand, < j < i — 1) is invertible in C[X]m because ^ when 
j < i. Then If* G (^)C[X] Mp - Thus 

(Si)C[X] Mp C (^)C[X] Mp . 

That completes the proof. 

Corollary 1 Lef notations be as above. For 2 < i < n, MZem(Si/ J(i — 1)) = 
MZero(S' J '/ J(i — 1)). 77zen we can rewrite eq.^2} as 

n 

MZero(S) = MZero(C/ J{n)) U |J MZero(5|/J(z - 1)). (3) 

i=2 

Now, we may compute, for each S[, a characteristic set and obtain a zero decom- 
position with multiplicity of MZero(S' l '). Finally, by the discussion before Proposition 
12 the process must terminate as long as = prem(/ l ei , C) = never occurs. Then we 
have a triangular decomposition with multiplicity of the original polynomial system. 

The main algorithm can be outlined below. 

Algorithm: ZeroDecompMulti 

Input: polynomial system PS which satisfies (PS) is a zero-dimensional ideal. 
Output: a zero decomposition with multiplicity of MZero (PS). 

1. Initialize SET1={ [PS, 1]}, SET2=0 and SET3=0. 

2. Estimate a bound on the number of zeros (counted with multiplicity) of PS, 
denoted by m. 

3. Repeat the following steps until SET1=0 

3.1 Choose [S, P] in SET1 and set SET1 = SET1 \{[5', P}}. 

3.2 For the polynomial system S, compute its characteristic set by Wu's method. 
If it's inconsistent, go to Step 3.1. Otherwise, denote the characteristic 
set by C = [C\(X), . . . , C n (X)], and let Ii be the leading coefficient of 

d(X), J(i) = n h (1 < * < n). 

k=l 
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3.3 Let SET2 = SET2 U{[C,Px J{n)}}. 

3.4 For alH (2 < i < n), compute n = premfP™, C). If r l = 0, let SET3 = 
SET3 U {[SUC, P x J(i ~ 1)]} and SET2 = SET2 \{[C, P x J(n)]}, else 
let SET1 = SET1 U{[SUCU r i; P x J{i - 1)]}. 

4. Output SET2 and SET3. 

If 7*j 7^ always occur, the rank of the characteristic set of new polynomial system 
must be less than that of the original one. So the termination of the algorithm is obvi- 
ous. Otherwise, we put the new polynomial system into SET3 as a part of the output. 
Eventually, we get a zero decomposition with multiplicity of PS: 

MZero(PS) = |J MZero(T/P) U |J MZero(5/Q) (4) 

[T,P]eSET2 [S,Q}GSET3 

Now we present an example to demonstrate the algorithm. 

Example 1 [3 ] Consider a system of 3 polynomials PS — {x 2 + y + z — 1, x + y 2 + 
z — 1, x + y + z 2 — 1}. We try to compute a zero decomposition with multiplicity of PS. 
First, we use the Bezout bound on the number of zeros, m = 8, which is the product of 
all total degree of polynomials in PS. Next, we compute a characteristic set of PS: 

C = [x 2 (x 2 + 2x - 1)0 - l) 2 , x 2 (x 2 + 2y - 1), x 2 (x 2 - 1 + 2z)\. 

According to Theorem 1, MZero(PS') = MZero(C/x 4 )UMZero(PS'l), where PSl 
PS U C U {a; 16 }- We continue to compute a zero decomposition ofMZero(PSl). The 
characteristic set of PSl is: 

CI = [x 2 , y - x - y 2 , 1 - z - y}. 

Therefore MZero(PS'l) = MZero(Cl). Finally, we get a zero decomposition with 
multiplicity of MZero( PS) as follows: 

MZero(PS') = MZero(C/a; 4 ) U MZero(Cl). 

Because rj = never happens during the process, the decomposition is triangular. 

Now, we know that PS can be decomposed with multiplicity into 3 branches: 

Ti = [x 2 + 2x - l,x 2 + 2y - l,x 2 - 1 + 2z], 
T 2 = \{x - l) 2 , x 2 + 2y - 1, x 2 - 1 + 2z], 
T3 = [x 2 ,y - x - y 2 , 1 - z - y]. 

The multiplicity of each zero in T\ is 1, and the zeros in Zero(T2) or Zero(T3) all have 
multiplicity 2. 

The following example points out that some polynomial systems cannot be decom- 
posed into triangular form through the algorithm ZeroDecompMulti. 

Example 2 [3 / Consider a system of three polynomials PS — {x 3 —yz^y^ — xZjZ 3 — 
xy}. The Bezout bound on the number of the zeros is m = 27. We compute a charac- 
teristic set of PS: 

C = [x 7 - x 11 ,yx 4 - yx 8 , -zx 5 + zx 9 }. 
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According to Theorem 1, we know MZeio(PS) = MZero(C/(x 4 — x s )(x 9 — x 5 )) U 
MZem(PS UCU(x 4 - x 8 ) 27 ) = MZero(C/(a; 4 - x a )(x 9 - x 5 ))UMZero(PS U C) 
since prem((a; 4 — x s ) 27 ,C) = 0. For PS1 = PS U C, its characteristic set is 
also C, so algorithm ZeroDecompMulti cannot get a further decomposition of 
MZero(PS'l). Hence, the algorithm gives us a decomposition ofWLZevo(PS) as fol- 
lows: 

MZero(PS) = MZero(C/(ir 4 - a; 8 )) U MZero(P5 U C). 

Actually, there is no zeros in MZero(C/(x 4 — a; 8 )), so the decomposition is mean- 
ingless. 

5 Improvements on the algorithm 

We discuss some possible improvements on the algorithm ZeroDecompMulti in 
this section. 

First, we consider the bound on the number of zeros of the original polynomial sys- 
tem, which is set to be a constant m in the algorithm ZeroDecompMulti. Obviously, 
the smaller the m is, the easier the computation will be. 

Intuitively, we can set m to be the number of zeros (counted with multiplicity) 
of the original polynomial system, which can be obtained by computing the Grobner 
basis of the system. But we will have some additional time cost for computing the 
Grobner basis, which violates our intention obviously. So obtaining a reasonable bound 
according to the form or characteristic of the polynomial system with little cost is more 
intelligent. 

Another strategy is to define m as a variable, and update it during the decompo- 
sition process. Review eq.©, suppose that the number of zeros of polynomial sys- 

n 

tem S does not exceed mo, then that of 1J MZero(5i/J(« — 1)) must be no more 

than mo — |MZero(C/J(n))|. Here |MZero(C/ J{n))\ means the number of zeros 
which satisfy C — and J(n) ^ (counted with multiplicity). As a result, we could 
make the following adjustment in the algorithm ZeroDecompMulti: when adding 
[PS,P] into SET2, set to = to - |MZero(P5/P)|. Usually, |MZero(P5/P)| is easy 
to calculate, so the computation will be more efficient. 

Second, we consider the situation that i\ = prem(i™ , C) = occurs. Since the 
characteristic set of computed by WuChar Set must be C, we could not get a further 
decomposition of MZeio(S' i / J(i — 1)) by the algorithm ZeroDecompMulti and 
only record this situation in SET3. We may try the following strategy. 

Assume that n = prem(/ l m , C) = and CjYj = CjjX) - h x x f c s(CdX)) ^ ^ 
then we replace Si in eq.(2) by Sf — S U C U {Ci(X) }. The proposition below 
proves that the multiplicities of zeros are still unchanged. 

Proposition 3 Let the notations be as above. Ifri = prem(/™\ C) — 0, then 

MZem(Si/J(i - 1)) = MZero(5 l tt / J(i - 1)). 
Proof. Since i\ = prem(Z™, C) = 0, it's easy to prove that 
MZero(S l /J(i - 1)) = MZero(5' UCU {IT}/J{i ~ 1)) = MZero(5' U C/J(i - 1)). 
Because Zero(5 U Cj J(i — 1)) = Zero^' / J(i — 1)), we know from Lemma[3]that 
MZero(5 U C/J(i - 1)) = MZero(5,| t / J(i - 1)). 
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That completes the proof. 

Now we can try to decompose Sf as long as Ci(X) ^ 0. Let 

7i=prem(apO"\C), 

it's easy to prove that the conclusion still hold if we replace Ci(X) in Proposition^ 
by r Then, the rank of the characteristic set of Sf must be less than that of C, as long 
as r'^ 0. 

Example 3 Consider a system of three polynomials 

S = {x 2 + y, Axy + 2y 2 , (x + y)z 2 + z + 1}. 

According to the algorithm ZeroDecompMulti, we compute a characteristic set of 
S: C — [— 2x 3 + x 4 , x 2 + y, (x — x 2 )z 2 + z + 1], and get a decomposition 

MZem(S) = MZero(C/x - x 2 ) U MZero(Si), 

where S x = S U C U {{x - x 2 ) 12 }. Since prem((x - x 2 ) 12 , C) = 512x 3 ^ 0, we 
compute a characteristic set of S[ — SL)G'tl{x 3 }: C\ — [x 3 ,x 2 +y, (x— x 2 )z 2 +z+l}. 

50 the decomposition of Si is 

MZero(Si) = MZero(Ci/x - x 2 ) U MZero(S , 2 ), 

where S 2 = S[ UCi U {(x-x 2 ) 12 }. Clearly, MZero(Ci/x - x 2 ) = 0. The decompo- 
sition process will terminate because prem((a: — x 2 ) 12 , C\) = 0. However, this is not 
a triangular decomposition. Now, we use the method based on Proposition^ replacing 

51 by S\ = S[ U Ci U {(z + l) 12 } before computing the characteristic set. Finally we 
get a characteristic set C 2 = [x 3 , x 2 + y, (-1 + 22x - 232x 2 )z - 1 + 21x - 211x 2 ] 
of S\ and a complete triangular decomposition of S: 

MZero(S') = MZero(C/x - x 2 ) U MZero(C 2 ). 

Finally, we take into account the case of redundant components in the zero de- 
composition. Obviously, the decomposition will have a redundant component if some 
leading coefficients of the characteristic set have the same prime factors. Therefore, 
we should factorize these leading coefficients and consider the factors one by one. The 
multiplicities of zeros remain unchanged as long as the condition in Lemma [3] holds. 
This strategy may let us get a better decomposition sometimes. 

Example 4 Review the polynomial system in Example 2: 

PS = {.t 3 — yz, y 3 — xz, z 3 — xy}. 

A characteristic set of PS is C — [a; 7 — x , yx 4 — yx s , — zx 5 +zx 9 ], the corresponding 
leading coefficients are I = [—1, x A — x 8 , —a; 5 + a; 9 ]. All the square-free factors in 
I are x,l — x, 1 + x, 1 + x 2 . So we know that PS has a zero decomposition with 
multiplicity as below: 

MZero(PS) = MZero(PS* U {x 27 })U 

MZem(PSU {(1 - a;) 27 }/a;)U 
MZero(PS'U {(1 + a;) 27 }/a;(l - x))U 
MZero(PS* U {(1 + x 2 ) 27 }/x(l - x 2 )). 
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After decomposing each part of it, we have the final result: 

MZero(PS') = MZero(PS' U {x 7 , x 4 y, x 5 z})LI 

MZero([l - x, -y 4 + 1, -y 3 + z])U 
MZero([a; + l,y 4 - 1, -y 3 - z])U 
MZero([l + x 2 , y 4 - l,y 3 - xz\). 

Although there is still one component which is not triangularized, we find 3 triangular 
sets that keep the multiplicities unchanged. This is a better result compared to that of 
Example 2. 
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